Absence of interaction corrections in graphene conductivity 
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Graphene [l[ has several peculiar properties originating 
from its perfect two dimensionality and from the Dirac- 
like nature of its charge carriers at half-filling. In partic- 
ular, recent optical measurements Q show that at half- 
filling and small temperatures, if the frequency is in a 
range between the temperature and the band-width, the 

conductivity is essentially constant and equal, up to a 

2 

few percent, to 00 = %- § ; such value only depends on the 
fundamental von Klitzing constant h/e 2 and not on the 
material parameters, like the Fermi velocity. These ex- 
perimental results confirm the theoretical predictions Q 
based on the description of graphene in terms of massless 
non-interacting Dirac particles [1, [j| ; lattice effects have 
been taken into account in Q . Since truly universal phe- 
nomena are quite rare in condensed matter (an example 
is provided by the Quantum Hall effect), it is important 
to understand whether this apparently universal value is 
just an artifact of the idealized description in terms of 
non-interacting fermions or rather it is a robust property 
still valid in the presence of electron-electron interactions, 
which are certainly present and expected to play a role in 
real graphene. This question is entirely analogous to the 
one concerning universality in the quantum Hall effect 
0, a notoriously difficult and still open problem. 

The effects of the electron-electron interactions on 
the graphene conductivity have been investigated in the 
Dirac approximation by perturbation theory both in the 
presence of long- and of short-ranged interactions; how- 
ever, lowest order explicit computations have produced 
different results [8l4l0j]. depending on the regularization 
scheme (momentum cut-off or dimensional regulariza- 
tion) chosen to cure the spurious ultraviolet divergences 
introduced by the Dirac approximation. In Q it was pre- 
dicted that in the presence of electrostatic interactions 
the zero frequency conductivity tends to zero, while in 
8, 10] it was argued that it converges to the free Dirac 
one, as a consequence of the divergence of the Fermi ve- 
locity; however, if screening or retardation effects are 
taken into account, the Fermi velocity is known to satu- 
rate at low frequency llHl3l|. in which case it is unclear 
what to expect. The extreme sensitivity of the conduc- 
tivity computation to approximations or regularizations 
(see also [l3|) calls for a rigorous analysis. 

In this paper we consider the Hubbard model on the 



honeycomb lattice, as a model of monolayer graphene 
with screened interactions. While in general the under- 
standing of the low temperature behavior of the Hubbard 
model is a formidable challenge for theoreticians, in the 
case of the honeycomb lattice at half filling the methods 
introduced in [12j and based on constructive Renormal- 
ization Group have proved to be quite effective. Using 
these techniques, we rigorously establish the exact (non- 
perturbative) vanishing of the interaction corrections to 
the conducivity in the zero frequency limit. All Feyn- 
man graphs contributing to the conductivity cancel out 
exactly in the limit, a statement analogous to the Adler- 
Bardeen theorem in quantum electrodynamics 0, EH • 

We introduce creation and annihilation fermionic oper- 
ators ipi = (ai,bt ? ) = £~ 2 Er cR e ±lfe for 
electrons with spin index a =J \^ sitting at the sites of 
the two triangular sublattices A^ and of a periodic 
honeycomb lattice of side L; we assume that A^ = A has 
basis vectors h t 2 = |(3,±v3) and that Ab — Aa + Sj, 
with £1 = (1,0) and 62,3 = h(— 1>±\/3) the nearest 
neighbor vectors; Ba = {k — n%Gi/L + n^G^jL ■ 
< Ui < L} with Gi, 2 = ^(1,±V3) is the first 
Brillouin zone (note that in the thermodynamic limit 

L- 2 E- keBA "> \B\- l f B &> with \ B \ = 8vr 2 /(3V3)). 
The grand-canonical Hamiltonian at half-filling is Ha = 
Hjl + UV\, where Ho is the free Hamiltonian, describing 
nearest neighbor hopping (t is the hopping parameter): 
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and Va is the local Hubbard interaction: 



x£A B cr=tl 



x£Aa T=ti 



(2) 

In order to define the lattice current and the conductiv- 
ity, we modify the hopping parameter along the bond 
(x, x + Sj) as t — > tg j (A) — t expjie A(x + sSj) ■ 5j ds}, 
where A(x) G M 2 is a periodic continuum field on 
S A = {x = L^ih + L&h ■ & G [0,1)}; its Fourier 
transform is defined as A(x) = | <S a | — 1 X^ 6 x> A ^-p e_lpa: , 



where |«Sa| 



3%/3 



L? and V A = {p = nxGt/L 
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/L : rii £ Z} (note that in the thermodynamic 
limit I^SaT 1 Y,fev A ~> ( 2?r )~ 2 Ir 2 d P)- If we denote by 
H(A) = Hl({t Sd (A)}) + UVa the modified Hamiltonian 
with the new hopping parameters, the lattice current is 
defined as jg — — , which gives, at first order in A, 



Jp — 3p + 



(2tt) 2 ^« 9 ' 



(3) 



where J is a shorthand for |<Sa| 1 J2qev A an( lj if 

J B 4fe is a shorthand for X^Ba anc ^ 7? 



J _ l-e" 



ipSj 



3p 



a,3 



dk 



- ie tY I ^( + 6 - ^-^-^) 



-b+ aZ ^nic+*(^^- ? 0) 
is the paramagnetic current and 



(4) 
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(5) 



is the diamagnetic tensor. The conductivity, at Matsub- 
ara frequency Pq € 27r/3 _1 (Z + i) and in units such that 
ft = 1 , is defined via Kubo formula as [6| 



a lm (POJ " K | 

Po|^a| 



(6) 



where, if 3 = Tr{e"^ KA }, (■> = 3- x Tr{e-^ A -} and 
Xo = e HAX °Oe~ HAXa for a generic operator O, 

Jo 

(7) 

It is known that in general the interaction modifies the 
values of the physical quantities; for instance, the Fermi 
velocity vf, the wave function renormalization Z and the 
vertex functions are known to depend explicitly on the 
interaction moreover, it was proven in [l2| that vp, 
Z and the vertex functions are analytic functions of U for 
\U\ small enough, uniformly as /3, |A| — > oo. In this Letter 
we prove a similar result for the conductivity. Moreover, 
we prove that in the thermodynamic, zero temperature 
and zero frequency limit, the conductivity is universal, 
i.e., it is exactly independent of U. 

Theorem. There exists a constant Uq > such that, 
for \U\ < Uq and any fixed po, erf' (po) is analytic in U 
uniformly in j3,A as /?, |A| — > oo. Moreover, 



liin + i™ J} m u lm^) = To Sk 

P0->0+ P—KXl |A|-»oo n I 



(8) 



Note that the limit j3 — > oo is taken before the limit 
Po — > + . In other words, the theorem says that the 



interaction corrections to the conductivity are negligible 
at frequencies /3 _1 <C po <C t. 

Proof. The idea of the proof is based on the two main 
ingredients: (i) exact lattice Ward Identities (WI) re- 
lating the current-current, vertex and 2-point functions; 
(ii) the fact that the interaction-dependent corrections 
to the Fourier transform of the current-current correla- 
tions are differentiable with continuous derivative (in con- 
trast, the free part is continuous and not differentiable 
at zero frequency). This last property follows from the 
non-perturbative estimates found in [12j, which we now 
briefly recall. The generating functional for correlations 
can be written in terms of a Grassmann integral: 



W(A,\) = / p( rfl A) e V(i/0+(V>,A)+B(A,V) 



(9) 



where, if k = (ko,k) with ko the Matsubara frequency, 
P(dtp) is the fermionic gaussian integration for ip^a ~ 
( a kff^kir)' with inverse propagator 



7 - x (k) = -Z 



iko vofl*(k) 
voQ(k) iko 



(10) 
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with Zo = 1, Vq = |i and f2(fc) = |X)j=i 2 : 
(note that <?(k) is singular only at the Fermi points k = 
k± - (0,3p,±^)). Moreover, B(A,4>) = 



E 
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(11) 



and (tp, A) = L dxo YlgeA^i + V'x ]• The response 
function K@' A (p) corresponds to the spatial components 
of the tensor K^(p) = 5A ^f 5Av{ _ v) W(A, 0)| A=Q , with 
/i, v = 0,1,2. Performing the phase transformation 
ip± -> e ±lea *ip± in Eq.@, we find 

W(A + da,Xe lea ) = W(A,X) , (12) 

which implies the following lattice Ward Identity [l6| 
2 



E?v^(p) = ° » 

ju=0 



(13) 



for all v £ {0, 1,2}. On the other hand, the functional 
integral Eq.® can be evaluated in terms of an exact 
Renormalization Group (RG) analysis, described in full 



detail in [12J. We decompose the field ip as a sum of 



fields ip( k \ living on momentum scales |k — k^ 
with h < a scale label; the iterative integration of the 
fields on scales h < h' < leads to an effective theory 
similar to Eq. ^ with a cut-off around the Fermi points of 
width 2 h and with a scale dependent propagator gr'-^(k) 



3 



with the same singularity structure as Eq. (fTO|) . with Z 
and vq replaced by Zh and Vh, respectively (the effec- 
tive wave function renormalization and Fermi velocity 
on scale h). Moreover, setting for simplicity A = 0, at 
scale h the interaction V(ip) + B(A,ip) is replaced by an 
effective interaction V(<Q (ip(<V) + B^(A, ip^ h )), with 
the effective potential y(-' 1 ) (ip(- h >) a sum of monomials 
in i])(- h > of arbitrary order, characterized at order n by 
kernels W^o(xi, ...,x„) that are analytic in U and de- 
cay super-polynomially in the relative distances |x; — Xj| 
on scale 2 ; moreover the effective source is given by 



2 

IM=0 



dp 



A M (p)j M (p) + BW (14) 



where j (p) = ~^J2a I J^W^k+p.^oik,?)^, 
Sij exp{-ipiS i2 } and 



T(k,p) 



3^ e +i(k+MSi~Si) 




(15) 



Finally, B^ is a sum of monomials in ( A, tp) of arbi- 
trary order, characterized at order n in tp and to in A by 
kernels W^%{x.i,. . . , x n ; yi , . . . , y m ) that are analytic in 
U, decay super-polynomially in the relative distances on 
scale 2~ h and are non-zero only if to > l,n > and 
to + n > 3; in particular, for all < 6 < 1, they satisfy 
the bounds (proved in [12[), 



J dx.2 ■ ■ ■ c£x„dyi . . . dy m \ wjfy (xi , . . . , x„; yi, . . . , y m )| 

< (const.)|e| m 2( 3 -"- m )' l (l - 6 m>0 + \U\2 6h ) . (16) 

The bounds Eq. (|16l) are non-perturbative (i.e., they are 
based on the convergence of the expansion for the kernels 
W^). They are obtained by exploiting the anticommu- 
tativity properties of the Grassmann variables, via a de- 
terminant expansion and the use of the Gram-Hadamard 
inequality for determinants, see 12j. The factor 2 0h in 
the bound will play a crucial role in the following and 
reflects the fact that the scaling dimension 3 — n — m is 
always negative for n > 2. The running coupling con- 
stants Zh,Vh, Zu,h satisfy recursive equations (beta func- 



related by Ward Identities. Indeed, proceeding as in (13| . 
we consider a reference model defined in a way similar 
to Eq.©, with the important difference that the Grass- 
mann integration P(dip) is modified into P>/ l (ciV'), whose 
propagator differs from the original one by the presence 
of a smooth infrared cutoff selecting the momenta > 2 h ; 
performing the phase transformation ip^ — > e ±la *ipx m 
this functional integral, we find the analogue of Eq. (fl"2|) . 
which implies 



Zo.h 



= l + 0(U2 8h ) , 



Z- 



2,h 



= 1 + 0{U2 6h ) 



Z h v h Z h v h 

where the corrections 0(U2 6h ) come from the symme- 
try breaking terms due to the infrared cut-off function. 
Therefore, the effective parameters are related by ex- 
act identities] the vertex density renormalization Z^ is 
equal, up to negligible terms, to the wave function renor- 
malization, and the current renormalization zj^' is equal 
to the product of the effective velocity and the wave func- 
tion renormalization 11711. 



We can write K^ v = K^' +K]X' , where the two terms 
in the right hand side correspond to the paramagnetic 
and diamagnetic contributions to K^, see Eq.©. Note 

that Kj£'(po,0) is independent of po; using Eq. (fT6|) . we 
find that [K^ipoM < (const.)|e| 2 E°=-oo 2 ' S . which 
is finite. On the other hand Kp, v (x) = 



(D) 



h=- 



(Z, 



lh) 2 J 



(2ny\B\ 



e ipx F ft (k,p). (18) 



•Tr{r M (k,fic h (k)r v (k+p, -f)C h (k + P )} + eW(x 



where the first term corresponds to the zero-th order in U 
in renormalized perturbation theory (Fh (k, p) is a suit- 
able smooth cutoff function constraining |k — k^| and 
|k-k±+p| to be~ 2 h and such that J2h=-oo F h{K P) = 
1; moreover, Z/ l C^ 1 (k) is given by Eq.(fT0f with Z a ,v a 
replaced by Zh,Vh) and, for all N > and suitable con- 
stants Cat, 



\H$( X )\<C N \U\ 



2(4+e)ft 
l + (2 h |x|)" 



(19) 



tion equations) that, due to the bound Eq.JIB}, lead to As compare d to the zero-th order contribution to K 



bounded and controlled flows, i.e., Z(U) — lim/ l _ i ._ 00 Z, 
Zn(U) = hmft^-oo Z^h and v F (U) = lim^-oo Vh are 
analytic functions of U , analytically close to their unper- 
turbed values Zq = Zo.o = 1 anf l Zi.o = ^2,0 = i'o = §*j 
see [12J. The analyticity of the kernels of the effective 
potential and of the h — > — 00 limits of the running cou- 
pling constants implies the analyticity of the imaginary- 
time correlation functions (see [HI]) and, similarly, the 
analyticity of c^' A (po) claimed in the main theorem. 

We are left with proving the universality result Eq. © . 
To this aim, it is important to notice that Zh,Vh, Z^j t are 



the dimensional bound on HfrJ has an extra factor 2 eh , 



following again from Eq. dTBl . From Eq. fT^) . 
I4^(x-y)|< (const.) ' 



1 



(20) 



that is, -fT^j^x) is absolutely integrable and, therefore, its 
Fourier transform in the thermodynamic and zero tem- 
perature limit is continuous at p = 0. Combining this 
remark with the WI Eq. fl3|) . we find that K^O) = 
0. In fact, setting, e.g., p 2 = 0, K n (po,px, 0) = 



4 



( — Po/pi)Koi(po,pi, 0); taking first the limit po — > and 
then pi — > in the right hand side, we get Kn(0) = 0; 
proceeding analogously, we find that K^ u (0) = for all 
/j, v e {0,1,2}. 



On the other hand Eq.(Q7l) implies that = v F (U) + 

O{U2 0h ), so that k£\x) = /^ ;0) (x)+^ ;1) (x) where 
K^'°\x) is the paramagnetic response function for the 
model with Hamiltonian H K {\v F {U)) and |^ ;1) (x)| < 

with < 9 < 1. There- 



(const.)|C/|(l + Ixl^ 8 )- 1 , 

(p-i) 

fore, the Fourier transform of iv^„' is differentiable 
in p and its derivative is continuous at p = 0. A 
similar decomposition can be performed in the diamag- 
netic term, so that, defining K^u 



and using the Wl Eq.®, (P 



from 



which, setting, e.g., p2 = 0, we obtain k[\\pq,pi, 0) = 
(j>o/pi)' 2 Kqq\po,pi, 0); deriving with respect to po both 
sides and taking first the limit po — ¥ and next p\ — > 
in the right hand side, we get 8 Po K[^(0) — 0; proceed- 
ing analogously, we find that d Pp K${0) — for all 
p,fj-,v G {0,1,2}. Note the crucial role played by the 
continuity of the derivatives of Kj^J , which allowed us to 
exchange the zero frequency and zero momentum limits, 
as compared to the order in Eq.©. 

In order to compute the conductivity, we are left with 
the contribution associated to a free theory with Fermi 
velocity equal to vf(U) that, for the 11 component, set- 
ting v = vf{U), reads: 



CTll 



lim lim lim erf q ' A (iOo) 

p ->-0+ /3->oo |A|->-oo 



dk^ f S(k)-S(k- 



2eV 



lim 

po->o+ 



dk Q 
~2~7 



Po 



L (fc,o)^(k)r 1 (fc,o)} 



The latter integral is not uniformly convergent in po] in 
particular, it is well known that one cannot exchange the 
limit with the integral [l4[ . The integral can be evaluated 
explicitly (using residues to compute the integral over ho) 
and leads to Eq. ([T3|) . A similar computation shows that 
C22 = en and that the off-diagonal terms vanish. ■ 

The above analysis can be extended to the case of 
long range electromagnetic interactions; in such case 
the wave function, density and current renormalizations 
have a strong (anomalous) power law dependence on the 
momentum and the Fermi velocity increases up to the 
speed of light [H EDI- Wis similar to Eq.(TT7]) are still 



valid and imply that, even if the effective parameters are 
strongly momentum dependent, the conductivity only de- 
pends weakly on the frequency in the optical range. 

In conclusion, we rigorously proved the non existence 
of corrections to the zero temperature and zero frequency 
limit of the graphene conductivity due to weak short 
range interactions. The proof is based on a combination 
of constructive Renormalization Group methods with ex- 
act lattice Ward identities. Remarkably, this is one of the 
very few examples of universality in condensed matter 
that can be established on firm mathematical grounds. 
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